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INTRODUCTION 


The  objective  of  this  program  is  to  apply  a  solid  state,  staring,  ultrasonic  hybrid  array  to  demon¬ 
strate  the  first  real-time,  portable  C-scan  diagnostic  ultrasormd  medical  imaging  system  with  three- 
dimensional  (3-D)  capability  for  battlefield  and  trauma  care. 

Current  B-scan  ultrasound  imaging  limits  the  visualization  of  3-D  structures  because  1-D  B- 
scan  arrays  only  image  2-D  slices.  These  slices  must  then  be  reconstructed  in  the  operator’s  mind  to 
obtain  a  mental  3-D  “picture”  of  the  item  of  interest.  Automatic  movement  of  the  transducer  with  com¬ 
puter  reconstruction  (offline)  has  recently  been  introduced  by  some  manufacturers,  but  the  systems 
are  slower,  more  costly,  and  harder  to  use  than  conventional  B-scan  imaging  instruments.  They  are 
not  well  suited  to  battlefield  use. 

Recognizing  the  importance  of  3-D  imaging,  the  Tiltrasound  industry  is  in  search  of  a  techmcal 
solution  to  building  high  resolution,  2-D  ultrasound  arrays.  Our  solution  is  unique.  It  consists  of  con¬ 
structing  large  imaging  ultrasoiuid  arrays  using  processes  derived  from  ARPA  supported  technology 
for  infrared  focal  plane  arrays  (IRFPAs).  We  have  directly  applied  existing  designs  and  processes  from 
thermal  imaging  systems  to  prove  our  ultrasoxmd  imaging  capability. 

This  technology  bridges  the  large  gap  between  taking  vital  signs  and  making  basic  observations 
in  the  field,  and  taking  CT  scans  back  in  a  hospital.  With  C-scan  technology  brought  to  the  field,  deep 
tissue  damage  assessment  is  done  quickly,  enabling  earlier  intervention.  This  is  true  for  either  battle¬ 
field  environments  or  EMT/ambulance  applications.  Real-time  imagery  also  affords  the  opportunity 
for  transmission,  remote  diagnosis,  and  archiving  of  information.  This  can  also  benefit  health  care  in 
remote  locations  where  specialists  and  expensive  imaging  diagnostic  equipment  is  not  available. 

.  Earlier  efforts  had  already  proven  the  feasibility  of  tibe  2-D  transducer  hybrid  array.  They  had  suc¬ 
cessfully  imaged  the  direct  illumination  of  a  focused  source  with  our  transducer  array  technology.  We 
believe  this  44  x  64  element  ultrasoimd  array  to  be  the  largest  operating  two  dimensional  (2-D)  ultra¬ 
sound  array  in  existence. 

During  Phase  I  of  this  project,  our  goal  was  to  demonstrate  real-time  C-scan  imagery  from 
reflecting  targets.  Our  approach  was  to  build  2-D  ultrasonic  hybrid  transducers  using  a  diced  array  of 
piezoelectric  elements  mated  to  an  existing  readout  integrated  circuit  (ROIC).  We  also  modified 
existing  hardware  to  enable  the  capture  and  analysis  of  real  time  C-scem  imagery.  lb  accomplish  this, 
we  modified  an  existing  IR  camera  to  provide  range  gate  and  timing  controls  necessary  for 
ultrasonography.  While  this  proved  to  be  an  expeditious  and  cost  effective  course  for  the  feasibility 
demonstration,  it  also  had  inherent  limitations. 

In  parallel  with  this  effort  and  in  support  of  the  real  time  imagery  Phase  I  goal,  the  IRFEA.-FM  pro¬ 
gram  estabhshed  a  transducer  fabrication  capabilify  based  on  the  IR  read  out  integrated  circuit  and 
infrared  focal  plane  fabrication  technology  developed  from  our  IR  detector  fabrication  technology  base. 

This  report  describes  the  work  performed  in  Phase  I  with  an  emphasis  on  the  key  technologies 
demonstrated,  together  with  the  challenges  overcome  and  those  identified  for  Phase  II.  As  described 
below,  the  Phase  I  effort  demonstrated  the  feasibility  of  an  ultrasonic  battlefield  imager  using  2-D 
hybrids  as  the  ultrasonic  transducer  by  imaging  reflected  energy  from  both  specular  and  diffuse  scat¬ 
tering  targets. 

Our  Phase  II  program  builds  on  the  success  of  the  initial  feasibility  components  developed  in 
Phase  I.  We  believe  our  approach  will  have  benefit  in  locating  foreign  object  fragments;  the  location 
and  extent  of  hemorrhage;  and  in  improving  techniques  for  guided  surgery  and  biopsy. 
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SCIENTIFIC  &  TECHNICAL  MERIT 


The  goals  of  our  medical  ultrasound  imaging  system  development  are  visualization  of  the  inte¬ 
rior  of  a  human  subject  with; 

•  Excellent  image  quality 

•  Real-time  frame  rates 

•  Image  manipulation  to  obtain  any  derived  view  or  cross-section 

•  Portability 

•  Low  cost 

Demonstration  of  imaging  with  a  2816  element  ultrasonic  hybrid  array  is  a  milestone  achieve¬ 
ment  in  medical  ultrasoimd. 

RELEVANCE  TO  AREA 

Our  imaging  ultrasound  array  is  the  application  of  dual  use  technology  derived  from  ARPA 
sponsored  developments.  Arrays  are  derived  from  IRFPA  technology,  using  acoustically  sensitive 
materials  in  place  of  the  IR  detectors.  Under  the  recently  awarded  flexible  manufacturing  program 
with  ARPA/MTO,  Loral  has  also  demonstrated  fabrication  of  ultrasound  hybrid  arrays.  By  using  the 
same  IR  hybrid  fabrication  line  that  supplies  DoD  and  commercial  IR  imaging  products,  ultrasound 
arrays  have  been  made  available  to  this  program  without  significant  additional  startup  expenses, 
providing  maximum  synergy  of  ARPA  activities. 

CAPABILITIES  AND  COMMERCIALIZATION 

We  have  assembled  a  team  that  includes  key  leaders  in  diagnostic  ultrasoimd  medicine,  equip¬ 
ment  development,  and  imaging  systems: 

•  University  of  Rochester,  Center  for  Biomedical  Ultrasound:  Dr.  Kevin  Parker,  Direc¬ 
tor.  The  center  is  the  oldest  and  largest  interdisciplinary  research  organization  in  the 
world  devoted  to  aU  aspects  of  medicsd  imaging. 

•  Loral  Infrared  &  Imaging  Systems  (LIRIS):  A  leader  in  military  and  commercial 
application  of  electro-optic  imaging  technologies,  LIRIS  will  design  and  fabricate  the 
ultrasonic  imaging  arrays  and  associated  electronics. 

•  Marvin  E.  Lasser,  Inc. :  Pioneer  of  3D  ultrasound  arrays  and  related  signal  and  image 
processing. 
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THE  BATTLEFIELD  ULTRASONIC  DIAGNOSTIC  IMAGER  (BUDI) 
FEASIBILnY  DEMONSTRATION  SYSTEM  DESCRIPTION 

The  Phase  I  feasibility  demonstration  system  is  a  bistatic,  staring,  ultrasonic  imager.  It  consists 
of  a  5  ]^z  transmit  transducer,  a  “physical  beamformer”  or  acoustic  imaging  mirror/lens,  a  two- 
dimensional  ultrasonic  hybrid  transducer  array,  and  camera  support  electronics.  These  elements  are 
described  below. 

The  2-D  Piezoelectric/ROIC  Hybrid  Array 

Figure  1  is  a  schematic  cross-section  of  the  hybrid  array.  It  consists  of  a  piezoelectric  plate  that 
has  been  partially  cut  (dicedXreticulated)  into  individual  array  elements  and  then  bonded  (hybrid¬ 
ized)  electrically  and  mechanically  to  a  silicon  integrated  circuit,  known  as  a  Read  Out  Integrated  Cir¬ 
cuit  (ROIC),  with  indixim  “bumps”. 


All  the  ROICs  used  in  Phase  I  were  designed  for  an  Infrared  Camera 
that  was  fabricated  eight  years  ago. 

This  existing  ROIC  has  136  x  192  individual  elements  on  50  fim  centers.  The  5  MHz  (center  fre¬ 
quency)  piezoelectric  array  was  designed  to  have  44  x  64  elements  on  150  pm  centers  so  that  each 
ultrasonic  element  is  connected  in  parallel  to  a  3  x  3  subarray  of  the  ROIC.  (Figure  1  is  a  concept  draw¬ 
ing  and  does  not  show  this  3x3  connection).  Since  the  acoustic  wavelength  in  water  at  5  MHz  is  300 
pm,  the  square  array  elements  are  thus  X/2  on  a  side. 

The  total  of  2816  active  array  elements  is  more  than  a  factor  of  20  larger  than  most  high  perfor¬ 
mance  commercial  ultrasonic  systems. 

ROIC  Operation 


A  schematic  block  diagram  of  the  ROIC  is  shown  in  Figure  2  together  with  the  horizontal  and 
vertical  readout  logic  which  allows  an  individual  array  element  (Figure  3)  to  be  addressed. 

Figure  4  is  a  block  diagram  of  the  entire  readout  system  and  Figure  5  describes  the  operation. 
Image  collection  begins  with  the  emission  of  an  ultrasound  pulse  from  the  transmitter.  A  negative 
range  gate  pulse  is  applied  to  all  2816  elements  in  parallel  through  the  detector  common  electrode  on 
the  front  surface  of  the  piezoelectric  material.  Its  delay  is  selected  to  coincide  with  ultrasonic  energy 
returning  from  the  object  plane  of  interest.  Its  width  is  also  rmder  user  control.  For  the  duration  of  this 
pulse,  a  parasitic  diode  in  the  amplifier  at  each  cell  acts  as  a  peak  detector,  storing  an  entire  frame  of 
ultrasoimd  information.  This  peak  detected  value  is  read  out  element  by  element  into  an  external 
transimpedance  amplifier  after  which  it  is  digitized  to  8  bits  as  Frame  A.  The  array  is  then  reset  by  the 
reset  pulse.  This  process  is  repeated  without  any  ultrasoimd  present  to  form  Frame  B,  which  is  stored 
in  memory  and  represents  the  “fixed  pattern”  noise,  due  to  different  electronic  offsets  or  other  noise, 
on  a  pixel  by  pixel  basis.  Frame  A  (ultrasound  present)  is  repeated,  however,  this  time  the  fixed  pat¬ 
tern  B  is  subtracted  from  the  values  of  frame  A  thus  improving  signal  to  noise  ratio.  The  (A-B)  frame 
is  then  formatted  and  output  to  a  gray  scale  monitor  as  conventional  RS-170  video. 
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Figure  1. 2D  Piezoelectric/ROIC  Hybrid  Array 
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Figure  2.  Ultrasound  Array  Readout  1C  (ROIC) 
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Figure  4.  Ultrasound  Camera  Video  Channel  Block  Diagram 
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Figure  5.  Experimental  Ultrasound  Camera  Response  and  Timing 

PHASE  I  PROGRESS  &  ACCOMPLISHMENTS 

PIEZOELECTRIC  ARRAY  FABRICATION 

The  piezoelectric  array  was  fabricated  from  a  wafer  of  PZT-5A  using  a  conventional  semiconduc¬ 
tor  dicing  saw  and  conventional  processes.  As  shown  in  Figure  6A,  the  first  arrays  were  "partially 
reticulated”,  i.e,  the  saw  cuts  were  penetrated  into  the  PZT  260  pm  (66%)  leaving  a  substantitd 
amoiint  of  material  at  the  front  (acoustic  port)  of  the  array.  The  saw  kerf  is  about  38  pm  and  the  aspect 
ratio  of  the  cuts  is  6.6:1.  In  these  first  arrays,  the  saw  kerf  was  left  unfilled.  Thirty  (30)  partially  retic¬ 
ulated  arrays  were  fabricated.  Experiments  disclosed  a  relatively  high  level  of  interelement  mechan¬ 
ical  crosstalk,  which  is  presumed  to  be  caused  by  the  partial  reticulation. 

Process  development  was  begun  to  produce  a  “fully  reticulated”  array  as  shown  in  Figure  6B. 
The  piezoelectric  material  available  for  this  work  had  a  two  layer  matching  section  bonded  to  the  front 
face.  Attempts  were  made  to  dice  completely  through  the  PZT  into  the  magnesium  layer  in  two  direc¬ 
tions,  however,  the  dicing  operation  caused  failure  in  the  PZT/magnesium  bond  resulting  in  loss  of 
array  elements  and  this  approach  was  abandoned. 

“Fully  reticulated”  arrays  were  produced  with  the  structure  shown  in  Figure  6C.  In  this  case  the 
PZT  was  sawed  in  one  direction,  back^led  with  epoxy,  then  sawed  at  right  angles  and  backfilled  a  sec¬ 
ond  time.  Finally,  a  silver  epoxy  electrode/matching  layer  was  applied.  Five  fully  reticulated  arrays 
have  been  fabricated.  Bowing  of  the  array,  presiamably  caused  by  shrinkage  of  the  silver  epoxy  layer 
has  been  a  significant  problem,  since  the  flatness  required  for  reliable  hybridization  is +5  pm.  Figure 
7  is  an  SEM  image  of  a  piezoelectric  array. 
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Figure  6.  Transducer  Array  Configurations 

Experimental  results  with  these  arrays,  although  limited,  have  confirmed  that  the  crosstalk  is 
substantially  reduced  with  the  fully  reticulated  structxu-e. 

HYBRIDIZATION 

The  hybridization  process  (Figure  8)  consists  of  metal  bumps  onto  the  back  side  of  the  array 
using  proprietary  IRFPA  techniques.  The  array  and  the  ROIC  (with  matching  indirnn  bumps)  are 
then  aligned  with  a  specialized  optical  probe  (Figure  8A),  the  probe  is  removed  (Figure  8B)  and  the 
two  parts  are  pressed  together  tmder  a  prescribed  force/time  profile.  There  is  no  heat  involved  in  the 
cold  welding  process. 

The  completed  hybrid  is  mormted  in  a  conventional  integrated  circuit  (PGA)  package  and  the 
ROIC  is  wire  bonded  to  the  pins  in  a  conventional  manner  (Figure  1). 
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CAMERA  DESCRIPnON 

As  part  of  the  work  in  Phase  I,  four  cameras  were  modified  to  add  the  range  gate  function  and 
timing  controls.  Additional  modifications  were  made  to  improve  dynamic  range  and  sensitivity. 

Details  of  camera  performance  are  shown  in  Table  1.  The  measxirements  were  made  with  the 
same  hybrid  array  and  eiq>erimental  arrangement  by  changing  the  drive  voltage  with  a  step-switched 
attenuator.  The  ultrasoimd  beam  was  refiected  f  rom  a  large  planar  reflector  which  was  imaged  by  the 
lens.  The  minimum  signal  levels  are  relative  numbers  for  intercomparison  of  the  cameras  only. 


Table  1.  Camera  Performance 


Camera 

1 

2 

3 

4 

Dynamic  Range  -  dB 

18 

28 

20 

28 

Minimum  Signal  -  dB 

-33 

-43 

-35 

-43 

Figure  7.  SEM  Image  of  Piezoelectric  Array 
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UPPER  CHUCK 


UPPER  CHUCK 


A 


B 


Figure  8.  Hybridization  Process 

IMAGING  EXPERIMENTS 

ParaUel  Approaches  to  Physical  Beamforming 

TWo  parallel  "physical  beamforming”  approaches  to  imaging  were  associated  with  the  Phase  I 
program:  an  acoustical  mirror-based  (reflective)  system  and  a  miilti-element  acoustic  lens  (refractive) 
system. 

A  schematic  drawing  of  the  mirror-based  system  is  shown  in  Figure  9.  A  spherical  cap  (virtual 
point  source)  transmitter  is  arranged  to  insonify  a  volume  such  as  the  body  through  a  membrane. 
Returning  echoes  &om  a  narrow  region  in  the  object  are  imaged  by  the  ultrasound  mirror  onto  a  piezo¬ 
electric  receiver  array,  which  is  controlled  by  the  ROIC.  After  processing,  the  acoustic  information  is 
converted  to  video  and  displayed  on  a  conventional  gray  scale  monitor. 

Figure  10  is  a  schematic  drawing  of  the  lens-based  system.  This  feasibility  system  is  immersed 
in  a  water  tank  and  also  uses  a  separate,  but  collimated  transmitter.  Returning  echoes  are  imaged 
with  a  two  element  acoustic  lens  onto  a  receiver  array  and  read  out  as  in  the  mirror  based  system. 
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Reflective  Imaging  System 

The  reflective  imaging  system  work  was  discontinued  after  initial  testing.  While  the  concept  has 
some  attractive  featxires,  this  particular  implementation  of  a  reflective  imager  suffered  from  a  num¬ 
ber  of  limitations: 


1.  Due  to  the  obscuration  of  the  receiver  array,  a  single  mirror  imaging  system  must  be  oper¬ 
ated  off-axis  which  substantially  increases  image  aberrations. 

2.  The  mirror  offers  only  a  single  optical  surface  for  correction  of  aberrations. 

3.  The  mirror  is  too  small  for  high  spatial  resolution.  It  should  be  fll.O  or  faster. 

4.  Diffraction  from  the  receiver  array  in  the  middle  of  the  ultrasonic  field,  the  mirror  and  sup¬ 
port  and  the  rubber  membrane  support  dominated  the  image. 

5.  Isolation  of  the  receiver  array  from  the  water  medium  was  attempted  through  a  latex  con¬ 
dom  which  presented  practical  difhculties  in  working  with  the  arrays  themselves  as  well  as 
an  uncontrolled  variable  in  the  acoustical  path. 


•  Imager  operational 

•  High  crosstalk 
observed  in  image 

•  Awaiting  fully 
reticulated  array  in 
LCC  package 


Figure  9.  Reflective  Imager  Setup 
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WATER  TANK  WALL 


Figure  10.  Refractive  imaging  Setup 


Refractive  Imaging  System 

The  refractive  imaging  system  was  designed  to  be  a  laboratory  feasibility  system  which  requires 
objects  to  be  immersed  in  a  water  tank  for  imaging. 

Transmitter 

The  transmitter  is  a  5  MHz,  25  mm  diameter,  unfocused,  commercial  off  -the-shelf,  piston  radi¬ 
ator.  The  tinfocused  sound  beam,  (shown  in  Figfure  11  as  a  CW  or  continuous  wave)  is  collimated  from 
the  face  of  the  transducer  to  approximately  500  mm.  In  actual  operation  the  transmitter  was  pulsed. 

Ultrasonic  Lens 

The  ultrasonic  doublet  lens  was  designed  using  standard  optical  lens  design  computer  software 
in  routine  use  at  LIRIS.  A  ray  trace  and  predicted  performance  parameters  of  the  uncoated  f'l.O  poly¬ 
styrene  lens  is  shown  in  Figure  12.  The  “as  designed”  85%  diffraction  blur  circle  diameter  numbers 
refer  to  object  space  and  also  include  aberrations.  Two  of  the  four  surfaces  of  the  lens  are  aspherical. 

The  lens  was  machined  using  a  numerically  controlled  lathe  to  a  nominal  tolerance  of  50  p.m  (or 
1/6  ultrasoimd  wavelength).  Measurement  of  the  lens  surfaces  with  a  depth  gage  accurate  to  ±  1  pm 
disclosed  that  three  of  the  four  surfaces  were  within  5  pm  of  the  design  value.  The  fourth  surface,  one 
of  the  aspherical  surfaces  proved  to  be  approximately  50  pm  from  the  design  value;  however,  it  was 
within  the  stated  machining  tolerance.  Tlie  “as  buUt”  numbers  in  Figure  12  reflect  this  change  in 
actual  lens  shape.  These  performance  differences  are  probably  unmeasurable  in  the  water  tank. 
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537  mm 
(Near  Reid) 


Far  Held 


179  mm 


25  mm  diameter 


Hgure  11.  IVansmIt  Beam  Profile  —  5  MHz-25  mm  Diameter 


Focal  Length  =  50  mm,  f/#  =  1.0 
Depth  of  Field  r  1  mm 
Minification  3:1 


FRACTIONAL 
OBJECT  HEIGHT 

85%  ENCIRCLED  DIFFRACTION  BLUR  (mm) 

AS  DESIGNED 

AS  BUILT 

0 

0.91 

0.92 

0.5 

0.87 

0.89 

1.0 

0.94 

0.96 

Figure  12.  Ultrasonic  Doublet  Lens 
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Shear  Waves  re  linage  Quality 

One  of  the  questions  posed  at  the  technical  reAdew  of  31  May  1995  was  the  influence  of  shear 
waves  formed  in  the  acoustic  lens  due  to  mode  conversion  of  the  compressional  (longitudinal)  at  the 
waterdens  interface.  Although  shear  waves  do  not  propagate  in  water,  they  can  be  generated  at  the 
water/lens  interface,  propagate  across  the  solid  material  of  the  lens  and  be  mode  converted  back  into 
compressional  waves  in  the  water  at  the  lens/water  interface.  This  spurious  wave  might  then  be 
detected  at  the  array  and  would  create  an  iinfocused  artifactual  signal,  limiting  contrast. 

Figure  13  is  a  ray  trace  of  such  an  image  using  the  shear  wave  index  of  refraction  in  polystyrene 
of  0.95.  The  shear  wave  velocity  of  sound  in  polystyrene  is  fortunately  very  close  to  that  of  the  sur¬ 
rounding  water  medium.  The  lens  thus  has  almost  no  “power”,  i.e.  rays  pass  through  essentially  unde¬ 
viated.  Another  way  to  express  this  is  that  the  lens  is  essentially  a  piece  of  “invisible”  material  as  far 
as  shear  waves  are  concerned. 

Since  the  lens  has  almost  no  power,  any  shear  waves  that  might  be  generated  are  spread  out 
across  the  lens  apertme.  From  the  ratio  of  the  areas,  the  array  thus  intercepts  at  most  about  1%  of  any 
shear  wave  energy  developed. 

A/R  Coatings 

Reflection  of  compressional  waves  at  the  lens/water  interface  is  a  more  serious  issue.  Figure  14 
is  a  graph  of  the  reflection  from  a  polystyrene  water  interface  as  a  function  of  incident  angle.  At  nor¬ 
mal  (0  degrees)  incidence,  the  reflection  coefficient  is  24.7%  at  each  interface,  which  is  rather  large. 
(Uncoated  optical  glass  typically  has  8%  reflection  at  an  air  interface.)  The  answer,  as  in  optical 
lenses,  is  anti-reflection  coating.  The  reflection  coefficient  of  a  simple  single  quarter  wave  coating  is 
also  shown  in  Figure  14.  It  can  be  seen  that  even  out  to  an  incident  angle  of  30  degrees,  this  simple 
layer  is  an  improvement  over  an  uncoated  lens.  This  simple  coating  has  less  than  1%  reflection  out  to 
about  10  degrees.  Figure  15  is  a  graph  of  the  frequency  response  of  such  a  coating  at  normal  incidence. 
A  2  MHz  bandwidth  can  be  achieved  rather  easily. 

Some  limited  initial  testing  was  done  with  one  A/R  coating  material  that  appeared  very  promis¬ 
ing.  It  was  not  successful,  however,  due  to  an  error  in  the  vendors  literature  for  the  elastic  modidus  of 
the  coating.  A/R  coatings  will  be  one  of  the  key  areas  for  development  in  Phase  II. 

Transducer  Hybrid  Interface 

Figure  16  is  a  schematic  drawing  of  the  hybrid  array  stage  with  the  moimted  PGA  package 
shown  partially  cut  away.  This  stage  mounts  the  hybrid  array  to  the  water  tank  as  shown  in  Figure  10. 
The  PGA  package  is  located  in  the  plastic  stage  on  four  edges  and  the  backside  and  is  held  in  place  by 
metal  clips,  one  of  which  is  shown  in  the  flgure.  The  wires  from  the  PGA  socket  are  led  out  through  the 
rear  of  the  stage. 

The  stage  itself  is  moimted  in  a  fixture  attached  to  and  penetrating  the  side  wall  of  the  water 
tank  (Figure  10).  A  stretched  membrane  moimted  on  one  side  of  this  fixture  contains  the  water  in  the 
tank  and  allows  the  array  to  be  acoustically  coupled  to  the  water  with  a  drop  of  oil.  This  arrangement 
has  proven  to  be  veiy  convenient  and  repeatable  for  testing  various  arrays. 

Display  and  Recording 

The  acoustical  image  is  displayed  on  a  black  and  white  monitor  without  any  additional  image 
processing  at  this  time.  A  3/4  inch  VCR  was  used  for  recording  images. 
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Amplitude  Reflection  Coefficient 


Shear  Wave  Index  of  Refraction  =  0.95 

Less  than  1  %  of  any  shear  wove  energy 
generated  Is  Incident  on  the  array 

Figure  13.  Ultrasonic  Doublet  Lens  —  Shear  Wave  “Image”  Is  out  of  Focus 


Angle  of  Incidence  -  degrees 


Figure  14.  Matching  Layer  •  Quarter  Wave  —  Water  to  Polystyrene  vs  Angle 


Amplitude  Reflection  Coefficient 


4,00  4.50  5.00  5.50  6.0i 

Frequency  -  MHz 


Figure  15.  Quarter  Wave  Matching  Layer  —  Water  to  Polystyrene/Normal  Incidence 


Figure  16.  Array  Stage  (Showing  Mounted  PGA  Package  Partially  Cut  Away) 
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Predicted  Performance  of  the  Refractive  Imager 

Figure  17  is  an  idealized  graph  describing  the  ultrasonic  signal  and  noise  levels  in  the  human 
body  as  a  function  of  target  depth  from  the  skin  surface.  It  is  based  on  the  average  value  of  attenuation 
in  the  body  of  1  dB/cm/MHz.  At  5  MHz,  the  attenuation  is  thus  5  dB/cm.  The  graph  is  referenced  to  0 
dB  which  is  the  reflection  from  a  perfect  reflector. 

A  perfect  reflector  at  the  skin  surface  thus  gives  a  signal  level  of  0  dB.  If  this  same  reflector  is 
located  at  a  target  depth  of  2  cm,  the  signal  level  is  reduced  by  20  dB  (4  cm  total  travel  distance  x  5  dB/ 
cm).  At  4  cm,  this  value  is  -40  dB,  etc. 

Targets  of  interest  in  the  body  are  generally  not  perfect  reflectors,  but  have  some  lower  reflec¬ 
tance.  Typically  large  specular  reflections  such  as  the  bovmdaries  between  fat  and  muscle  are  -20  dB 
re  a  perfect  reflector.  Non-specular  targets  such  as  blood  vessels,  tumors,  small  foreign  bodies,  tissue 
parench3mia,  blood,  fluid  filled  cysts,  etc.  range  in  scattering  power  from  nearly  in  the  specular  range 
to  extremely  small  values  depending  on  their  size  and  orientation.  In  the  absence  of  blood  vessels,  car¬ 
tilage  and  connecting  tissue,  a  large  relatively  uniform  tissue  such  as  liver  parench3mia  is  generally 
considered  to  be  a  “volumetric”  scatterer,  i.e.,  it  is  composed  of  a  statistical  ensemble  of  unresolved 
scatterers. 

Echo  levels  in  the  body  at  one  depth  cover  a  wide  amplitude  range,  for  example,  40  dB.  Any  imag¬ 
ing  system  attempting  to  image  from  the  skin  to  deep  inside  the  body  at  5  MHz,  thus  must  contend 
with  a  huge  d3mamic  range  on  the  order  of  100  dB. 

(assuming  average  value  of  1  dB/cm/MHz  attenuation) 


Target  Depth  -cm 

Figure  17.  Signal  and  Noise  Levels  in  the  Body 
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There  are  two  limits  to  the  ability  of  the  imaging  system  to  detect  smzdl  targets.  The  first  is  the 
electronic  noise  level  as  shown  in  Figure  17.  Clearly,  with  ordinary  techniques,  signals  lower  in  ampli¬ 
tude  than  this  noise  level  will  not  be  seen  at  any  depth.  The  electronics  noise  limit  is  generally  the  lim¬ 
iting  factor  at  a  certain  depth  (beyond  8  cm  in  Figure  17).  Targets  with  larger  scattering  power, 
however,  can  be  seen  out  to  greater  depths  (almost  14  cm  in  this  example). 

Increasing  the  input  acoustic  levels  can  boost  echo  signals  above  the  electronics  noise  level  and 
make  it  possible  to  image  deeper.  Practical  limits  in  the  transmitter,  presumed  safety  limits  and/or 
governmental  regulatory  limits  do  not  permit  acoustic  input  levels  to  the  humzm  body  to  be  raised 
beyond  a  certain  point. 

The  other  limit  is  the  target  clutter/contrast  ratio.  This  limit  is  determined  both  by  the  imaging 
system  itself  and  the  scattering  properties  of  the  tissue.  This  is  a  topic  of  great  contemporary  interest 
and  research  in  the  medical  ultrasound  field,  where  it  is  often  discussed  in  terms  of  aberration  correc¬ 
tion.  Unlike  the  case  of  the  electronic  noise  level  limit,  it  is  not  possible  to  increase  the  input  acoustic 
level  to  improve  target  contrast  ratio  because  it  is  essentially  constant  -  independent  of  input  levels. 
Note  that  the  target  contrast  ratio  applies  to  any  kind  of  beamforming,  including  “physical”  beam- 
formers  such  as  lenses  or  mirrors  as  well  as  electronic  beamforming. 

Reference  1,  entitled  “Scattering  from  a  Multiple  Random  Phase  Screen  Model  of  a  Random 
Inhomogeneous  Medium”,  by  Bernard  Steinberg  of  the  University  of  Pennsylvania,  is  a  reprint  of  a 
relevant  research  paper.  It  is  included  as  the  Attachment  to  this  report.  Quoting  from  this  reference: 

“Wavefronts  distort  when  ultrasound  propagates  through  an  inhomogeneous 

medium,  lowering  target  contrast  in  the  image  and  introducing  image  artifacts. 

Isotropic  scattering  reduces  the  contrast... 

“The  paper  presents  a  theory  of  the  effect  of  scattering  upon  important  image 

quality  parameters  such  as  target  to  clutter  power  ratio  (TCR)  and  lateral  resolu¬ 
tion.  Several  elements  of  the  theory  are  tested  against  experimental  data. 

“The  main  findings  are: 

•  When  target  strength  is  stronger  than  background  scattering,  the  principal 
scattering  mechanism  is  primary  (i.e.  simple,  single  or  nonmultiple)  scattering. 

•  Scattering  produces  a  halo  in  the  image  centered  on  each  scatter  source  or  tar¬ 
get.  Its  width  is  the  order  of  2-6  deg  at  3  MHz  and  (its  angular  extent)  is  propor¬ 
tional  to  wavelength,  (italicized  comment  added) 

•  TCR  within  the  halo  is  linear  with  2-D  array  area,  inversely  proportional  to  time 
delay  variance  per  unit  distance  and  propagation  length  and  decreases  as  the 
square  of  the  frequency  and  correlation  distance.  The  formula  is: 


TCR 


D(o^(ji-X  /D  ) 

^  ^  xm  m'' 


(Eq.  1) 


where  L  =  array  size 


D  =  propagation  path  length  or  tissue  thickness 


(0  =  27i{  is  the  angular  frequency 
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It-  Xom~  ^  ~  typical  scatterer  size 
^-nn  =  time  delay  variance 

Djn  =  tissue  thickness  in  which  is  measured 

“Target  dynamic  range  is  about  an  order  of  magnitude  smaller  than  TCR... 

“The  numerator  parameter  and  first  two  denominator  parameters  are  design 
parameters.  The  next  two  denominator  parameters  are  measured  properties  of 
the  tissue.  The  last  is  the  length  of  the  tissue  seimple  in  which  tissue  parameter 
measurements  are  made.  Thus,  time  delay  variance  and  correlation  distance  are 
the  primary  tissue  measimements  required  for  prediction  of  TCR  in  a  purely  scat¬ 
tering  medium. 

“The  lateral  resolution  was  found  to  remain  A/L,  the  width  of  the  diffraction  pat¬ 
tern,  zind  unaffected  by  the  scatter  field  so  long  as  the  scattered  energy  does  not 
obscure  the  main  image  lobe.” 

Refractive  Imager  Performance 

Images 

Figure  18  is  an  image  of  a  set  of  different  size  Allen  wrenches  in  a  plastic  holder  taken  with  one 
of  the  cameras  and  a  partially  reticulated  array.  The  wrench  sizes  and  spacing  between  them  are 
shown  in  Table  2. 
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Figure  18.  Ultrasonic  Image  of  Allen  Wrench  Set 
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Table  2.  Allen  Wrench  Set  Dimensions 


Nominal 

Size 

Width  of  Face 
(mm) 

Spacing 
Between  Faces 
(mm) 

9/64 

2.0 

3.9 

1/8 

1.8 

3.8 

7/64 

1.6 

3.5 

3/32 

1.3 

3.4 

5/64 

1.1 

3.2 

1/16 

0.9 

3.0 

0.05 

0.7 

The  wrenches  were  rotated  in  the  image  plane  for  msiximum  image  brightness  and  present  large 
specular  targets. 

The  field  of  view  is  approximately  25  mm  which  is  consistent  with  the  diameter  of  the  transmit¬ 
ter  transducer  beam  width.  The  wrenches  are  all  completely  and  easily  resolved.  Spatial  resolution  is 
thus  substantially  better  than  3  mm  and  as  scaled  from  the  photograph  is  estimated  to  be  about  1  mm, 
which  corresponds  to  the  theoretical  predictions  of  Figure  12. 

An  image  of  the  backbone  of  a  formalin  fixed  mouse,  made  at  the  University  of  Rochester,  is  shown 
in  Figure  19.  The  ossifications  of  the  vertebral  bodies  have  relatively  large  scattering  cross-sections. 
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Figure  19.  The  Backbone  of  a  Formalin  Fixed  Mouse 
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Quantitative  Array/Camera  Sensitivity  and  Minimum  Detectable  Signal 

Figure  20  shows  the  camera  sensitivity  and  minimxim  detectable  signal  (noise  level)  used  for  these 
images  superimposed  on  the  graph  of  Figure  17.  With  the  limited  sensitivity  and  dynamic  range  of  this 
camera,  it  is  clear  that  only  large  targets  close  to  the  surface  or  with  no  intervening  tissue  can  be  imaged. 

After  the  technical  review  of  25  September  1995,  camera  d3mamic  range  and  minimum  detect¬ 
able  signal  level  were  substantially  improved.  Figure  21  is  an  ultrasonic  image  of  parallel  nylon 
monofilament  line  targets  in  a  tissue  equivalent  test  object  (Model  539,  ATS  Laboratories,  Bridgeport, 
CT).  These  line  targets  are  0.12+5%  mm  in  diameter.  T^ey  are  spaced  in  5  mm  increments.  The  target 
in  the  center  of  the  image  in  Figure  21  is  7  mm  below  the  surface  of  the  test  object. 

Since  these  line  targets  are  less  than  X/2  in  diameter,  they  are  cylindrically  isotropic  scatterers. 
Their  image  is  thus  the  line  spread  function  of  the  image  system.  From  the  5  mm  target  spacing,  a  line 
width  of  approximately  1  mm  may  be  observed. 

While  it  is  gratifying  to  see  these  images  euid  the  recent  improvements,  several  points  must  be 
raised  regarding  system  sensitivity,  d3mamic  range,  and  electronic  noise  levels: 

1.  Array  sensitivity  is  relatively  poor.  This  is  due  to  many  factors: 

a.  The  tank  membrane  coupling  the  front  face  of  the  array  to  the  water  has  not  been  opti¬ 
mized. 

b.  The  acoustical  properties  of  the  array,  including  the  matching  layer(s),  mounting,  back¬ 
ing,  saw-cut  fill  material,  etc.  have  not  been  optimized. 

c.  The  indium  bump  mounting  of  the  ROIC  to  the  array  is  not  optimized. 

d.  The  ROIC  itself  was  designed  for  p3Toelectric  IR  detection,  not  ultrasound. 

2.  Camera  noise  levels  are  relatively  high 

a.  Electronic  noise  is  the  primary  limitation  to  detection  of  the  low  level  ultrasound  signals 
from  tissue  or  the  tissue  equivalent  test  object.  There  are  many  different  sources: 

Fixed  pattern  noise.  Array  elements  are  not  present  or  not  electrically  connected. 

Oscillations.  Low  frequency,  flicker  oscillations  in  the  image  are  a  cosmetic  defect. 

b.  Wbration  noise.  The  piezoelectric  array  has  response  at  low  frequencies  which  leads  to 
sensitivity  to  vibration  from  fans,  motors,  etc.  Use  of  a  vibration  isolated  optical  bench 
has  ameliorated  this  problem. 

c.  Artifacts.  There  are  various  artifacts  in  the  video  image. 

3.  Dynamic  range  is  still  limited  to  about  28  dB  best  case 

The  first  arrays  were  only  partially  reticulated  (partially  diced  through)  producing  relatively  high 
levels  of  cross  talk  in  the  array.  Later  arrays  with  complete  reticulation  appear  to  have  much  improved  iso¬ 
lation,  however,  they  suffered  from  other  defects  due  to  bowing  of  the  ceramic/silver  epoxy  matching  layer. 
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log  Signal  Level 

dB  re  Perfect  Reflector 


(assuming  theoretical  array  transducer  sensitivity) 


Target  Depth  -cm 

Figure  20.  Signal  and  Noise  Levels  in  the  Body  —  Phase  I  ROIC 


Figure  21.  Image  of  Monofilament  Nylon  in  Test  Object 
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Figure  22.  Signal  and  Noise  Levels  in  the  Body  —  Phase  II  Electronics 


Further  improvements  in  signal  to  noise  ratio  together  with  larger  ultrasonic  input  levels  will  be 
required  for  imaging  in  the  body.  This  will  be  achieved  by  designing  an  ROIC  specifically  for  the  ultra¬ 
sound  system  with  suitable  performance  (Figure  22).  Little  or  no  additional  effort  will  be  expended  to 
further  improve  the  existing  IR  ROIC  or  the  camera  electronics. 
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CONCLUSIONS 


Phase  I  Accomplishments 

Phase  I  of  this  effort  demonstrated  the  feasibility  of  an  ultrasonic  imager  emplo3dng  two-dimen¬ 
sional  hybrid  transducers.  Significant  demonstration  included  the  fabrication  of  two-dimensional 
piezoelectric  arrays  of  up  to  44x66  elements;  the  electrical  interconnection  of  these  arrays  to  sUicon 
readout  integrated  circuits;  the  fabrication  of  an  acoustic  lens;  and  the  imaging  of  specular  and  diffuse 
scattering  targets.  Images  collected  to  date  include  specular  reflections  off  metal  bar  targets,  diffuse 
reflections  off  monofilament  nylon  embedded  in  tissue  simulating  material,  and  the  backbone  of  a  for¬ 
malin  fixed  mouse.  The  system  collects  and  presents  these  images  in  real-time  (30  Hz  frame  rate),  and 
presently  resolves  targets  of  ~1  mm  at  a  distance  of 200  mm. 

OPPORTUNITIES  AND  CHALLENGES  IDENTIFIED 

Anti-Reflection  Coatings  for  the  Ultrasonic  Lens 

Anti-reflection  coatings  for  the  ultrasonic  lens  are  an  important  part  of  the  development  work  in 
Phase  II.  Polyst3n*ene  remains  the  material  of  choice  for  the  lenses  because  it  has  a  very  usable  index 
of  refraction,  relatively  low  ultrasonic  attenuation  and  an  acoustic  impedance  close  to  that  of  water.  It 
is  expected  that  only  a  single  quarter  wave  matching  layer  will  be  required  on  each  surface.  Selection, 
fabrication  and  performance  verification  of  this  layer  will  be  a  high  priority  for  the  early  part  of  Phase 
II  work. 

Electronics  Noise 

The  ROIC's  and  cameras  used  in  Phase  I  were  intended  for  a  different  application  and  compro¬ 
mised  the  result  that  were  achieved.  During  Phase  II,  these  components  will  be  specifically  rede¬ 
signed  to  exhibit  the  low  noise  levels  required  for  the  ultrasound  imaging  application  The  design, 
fabrication  and  integration  of  integrated  circuits  and  sensors  is  the  key  technology  strength  of  LIRIS. 

Target  Contrast  Ratio  Optimization 

In  the  Phase  I  final  review  of  25  September  1995,  the  issue  of  target  contrast  ratio  was  discussed. 
While  it  is  clear  that  this  is  a  limiting  factor  in  the  performance  of  any  ultrasonic  imaging  system, 
quantifying  this  parameter  for  an  ultrasonic  lens  based  system  is  another  critical  task  in  the  early 
part  of  Phase  II.  Methods  for  optimizing  this  parameter  such  as  apodization  and  other  sigru^l  process¬ 
ing  will  also  be  examined. 

The  results  of  this  work  will  be  compared  with  contemporary  electronic  beam  forming  methods. 
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ATTACHMENT 


“Scattering  from  a  Multiple  Random  Phase  Screen  Model  of  a 
Random  Inhomogeneous  Medium” 
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ABSTRACn* 

The  decrease  in  image  target  contrast  due  to  scattering  of  ultrasound  energy  by  an 
inhomogeneous  body  is  calculated  and  tested.  The  analysis  is  in  terms  of  measurable 
parameters  of  the  medium  (thickness,  phase  or  time-delay  variance  and  correlation  length)  and 
parameters  of  the  scaiming  equipment  (ftequency  and  array  size).  The  propagation  medium  is 
modeled  as  a  sequence  of  thin,  random  phase  screens  and  the  number  of  screens  is  related  to  the 
sample  thickness  and  to  the  size  of  the  scattering  bodies. 

It  is  found  that  (1)  scattering  produces  a  halo,  centered  on  each  target,  with  diameter  (2- 
6P  at  3  MHz)  proportional  to  wavelength,  (2)  target  contrast  within  the  halo,  defined  as  target  to 
clutter  power  ratio  (TCR),  increases  linearly  with  2-D  array  area,  is  inversely  proportional  to 
tissue  sample  thickness  and  time  delay  variance  per  unit  distance  in  the  tissue  and  decreases 
quadratically  with  correlation  distance  and  frequency  and  (3)  multiple  scattering  plays  no 
significant  role. 

Several  elempots  flfthe  thcoi^  are^tested  against  experimental  data.  Although  the 

amount  of  relevant  data  is  small,  the  theory  checks  well  against  what  does  exist. 


PACS  numbers;  43.80.Cs,  43.80£v,  43.80. Vj 

GLOSSARY  OF  TERMS 
L  array  size 

X,  f,  0)  =  2Ttf  wavelength,  frequency,  angular  frequency 

X/L  lateral  resolution 

phase  variance  per  screen 


*  Submitted  to  JASA.  Principally  supported  by  grants  from  NSF  (BCS92-09680)  and  the  Army  (DAMD17-93- 
J-3014). 
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time  delay  variance  per  screen 

xo 

e'l  correlation  distance  of  tissue 

Aus  =  X/tcxq 

half-power  width  of  the  scatter  spectrum: 

a  *  Jixo 

typical  scatterer  size 

TCR 

target  to  clutter  power  ratio  within  halo 

TDR 

target  dynamic  range  within  halo 

D 

propagation  path  length  or  tissue  thickness 

N-D/fl 

number  of  screens 

2  2 
^^41111,  Oxm 

measured  phase,  time  delay  variance  in  tissue 

2  2 

tissue  thickness  in  which  04)01  or  C^m  are  measured 

INTRODUCTION 

Wavefronts  distort  when  ultrasound  propagates  through  an  inhomogeneous  medium, 
lowering  target  contrast  in  the  image  and  introducing  image  artifacts**^.  Isotropic  scattering 
reduces  the  contrast  The  artifacts  are  discrete  multipath  arrivals  that  look  like  targets®.  The 
motivation  for  most  studies  of  this  topic  is  aberration  correction  for  large  aperture,  2-D  high 
resolution  imaging.  This  paper  deals  with  isotropic  scattering  only.  Discrete  multipath  will  be 
dealt  with  later.  The  paper  presents  a  theory  of  the  effect  of  scattering  upon  important  image 
quality  parameters  such  as  target  to  clutter  power  ratio  (TCR)  and  lateral  resoludon.  Several 
elements  of  the  theory  are  tested  against  experimental  data. 

The  main  findings  are: 

®  When  target  strength  is  stronger  than  background  scattering,  the  principal  scattering 
mechanism  is  primary  (i.e.,  simple,  single  or  nonmuldple)  scattering. 

®  Scattering  produces  a  halo  in  the  image  centered  on  each  scatter  source  or  target  It's  width 
is  the  order  of  2-6  deg  at  3  MHz  and  is  proportional  to  wavelength. 

o  TCR  within  the  halo  is  linear  with  2-D  array  area,  inversely  proportional  to  time  delay 

variance  per  unit  distance  and  propagation  path  length  and  decreases  as  the  square  of  the 

l2 

frequency  and  correlation  distance.  The  formula  is  TCR  «  - = - = - 

Da)2(7CXOin)^(tJTm/Uin) 

o  Target  dynamic  range  is  about  an  order  of  magnitude  smaller  than  TCR. 

§1  models  the  inhomogeneous  medium  as  a  cascade  of  thin,  random  phase  screens  and 
the  incident  radiation  from  the  target  as  a  plane  wave.  In  §2  the  scattered  power  Pscat  per  screen 
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is  calculated  as  a  function  of  the  phase  variance  per  screen,  and  also  the  ratio  of  Pjcat  to  the 
power  Ppw  of  the  emergent  plane  wave.  We  also  calculate  the  angular  width  of  the  scatter  Held. 
From  these  we  determine  the  TCR  in  the  target  image  and  how  the  lateral  and  contrast 
resolution  or  target  dynamic  range  are  affected  by  the  scatter.  The  target  is  the  source  of  the 
plane  wave  and  the  clutter  is  the  scatter  or  background  energy.  In  §3  we  calculate  the  effective 
number  of  screens  making  up  a  large  body  and  show  how  to  combine  the  multiple  scattering 
effects  between  the  many  screens.  We  then  calculate  the  quantities  above  for  the  entire 
propagation  medium.  §4  discusses  related  topics.  Summary  and  Observations  is  §5. 

1.  THE  MODEL 

When  a  point  target  or  source  is  imaged  through  a  weakly  scattering  medium  (small 
scattering  coefficient)  the  scattered  energy  creates  a  background  field  that  clusters  about  the 
target  image  and  decreases  its  contrast  with  the  background  and  therefore  its  visibility.  It  is  the 
same  phenomenon  that  produces  the  halo  around  stars  and  the  ring  around  the  moon  that  one 
often  sees  on  humid  nights.  The  ultrasound  halo  is  found  to  be  about  2-6  deg  wide  at  3  MHz 
centered  about  the  target.  Its  effect  is  minimal  with  small  transducers  whose  beams  are 
comparable  to  or  larger  than  the  halo.  But  with  narrow  beamwidth,  high  resolution  arrays  the 
scattered  energy  can  materially  reduce  target  contrast  The  calculations  use  only  measurable  or 
known  quantities.  The  scattered  energy  is  variously  called  scatter  or  clutter.  The  diminution  in 
TCR  in  the  image  is  calculated  as  well  as  the  angular  width  of  the  scatter  field  and  the  effects 
upon  lateral  and  contrast  resolution  (or  target  dynamic  range)  of  the  imaging  system.  The 
procedure  divides  the  medium  into  thin  slices,  calculates  the  relevant  statistical  properties  per 
slice  and  then  recombines  the  slices. 

A  unit  amplitude  plane  wave  from  an  impulsive  source  propagates  from  the  left  in  the 
positive  z  direction  through  a  random  inhomogeneous  medium,  defined  by  a  uniform,  isotropic 
distribution  of  scatterers,  located  in  the  x,z  plane  (Hgure  la)f  .  Each  slice  is  treated  as  a  thin, 
random  phase  screen  with  independent  scattering  occurring  in  each  screen.  Each  screen  is 
infinite  in  length  (x  axis).  Its  thickness  (z  axis)  is  the  order  of  the  scale  of  the  compressibility 
and  density  variations  of  the  medium,  denoted  a  in  the  paper.  Each  screen  is  defined  by  a  zero- 
mean,  Gaussian  random  phase  shift  process  <l)(x,z)  with  phase  variance  ^(x,z)  is  the  phase 

delay  at  frequency  f  within  the  signal  spectrum  through  a  thin  vertical  slice  perpendicular  to  the 
direction  of  propagation.  The  propagation  delay  t(x,z)  is  related  to  the  phase  delay  by  <>(x,z)  = 
(Dx(x,z),  CO  =  2jcf.  The  variation  in  t(x,z)  may  be  attributed  to  the  local  medium  peturbations,  as 
modeled  in  the  phase  screen  detail  of  Figure  lb.  The  scale  of  the  perturbations  is  denoted  a. 


■  Although  the  analysis  is  based  upon  a  plane  wave  source,  the  results  are  also  applicable  in  the  focal  plane  of  a 
focused  system  close  to  the  focal  point.  Sec  §4.1. 
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The  time  delay  variance  per  slice  cf  =  epofi.  The  correlation  function  of  the  phase  R(x)  = 
a|  p(x)  is  modeled  as  an  exponential  with  correlation  distance  x©.  Xq  is  related  to  scatterer  size 
a  by  fl  =  tcxq  (§2.5). 

Tissue  parameters  and  x©  are  obtained  from  measured  data.  Xq  is  the  value  of  x  for 
which  p  =  e'l.  Measurements  of  variations  of  arrival  time  across  wavefronts  disclose  ol. 
Measurements  are  usually  made  through  tissue  thickness  D  >  a,  number  of  slices  N  >  1.  The 
phase  variance  through  N  slices  is  a| and  the  relation  to  the  single-slice  variance  is  = 
No^  because  of  the  assumed  independence  of  the  phase  from  screen  to  screen  (§3.5).  The 
measured  value  of  x©  is  not  a  first-order  function  of  N. 

2.  SCATTERING  FROM  A  RANDOM  PHASE  SCREEN 

The  far-field  radiation  f(u)  from  a  single,  infinite  screen  (Figure  lb)  contains  an 
emergent  plane  wave  fpw(u)  plus  the  scattered  field  fscai(u).  The  plane  wave  is  reduced  in 
amplitude  to  account  for  the  energy  in  the  scattered  wave  as  it  passes  through  the  screen.  f(u) 
is  the  Fourier  transform  of  <)>(x),  in  units  of  the  reduced  angular  variable  u=sin9: 

f(u)  =  Jel<Kx)ejkxudx 

—  fpw(u)  +  fscat(u)  (I) 

k  =  1tU\  is  the  wavenumber  and  (|>(x)  is  the  random  phase  shift  previously  described.  f(u)  is 
decomposed  into  the  sum  of  a  plane  wave  fpw(u)  =  A5(u)  emerging  from  the  screen  and  a 
scattered  wave  fscat(ti)>  The  amplitude  of  the  plane  wave  is  A. 

2.1  The  emergent  plane  wave 

To  calculate  fpw(u),  f(u)  is  expanded  in  a  power  series  and  its  expectation  E[o]  taken. 

f(u)  =  J(1  +  j<|>  -  <J>2/2  -  j4>3/3!+  4,4/41 ..)  eilau  dx  (2) 

The  first  term  is  the  incident  plane  wave.  The  remainder  are  the  scattered  terms.  Their  mean 
values  (expectations,  spatial  dc)  add  (subtract)  coherently  to  (from)  the  incident  plane  wave  and 
the  total  establishes  the  amplitude  A  of  the  emergent  plane  wave.  Thus 

fpw(u)  =  Ef(u)  =  1(1  +  E[j(J)  -  4)2/2  -  j4>3/3!+  4)4/4!  •••])  ei*“u  dx  =  jAeJT«“  dx  (3) 

A  =  1  +  E[j(i)  -  4)2/2  -  j4)V3!+  4)^/4!  •••] 

=  1  +  E[i:n=l  jWn!] 

=  1  -  o2/2!  +3c^/4!  -  5o3a6/6!  +  7o5o3o|/8! 


(4) 


because  the  assumed  pdf  is  Gaussian  (all  odd  moments  are  zero). 

The  power  Ppw  in  the  emergent  plane  wave  is  proportional  to  A^.  Ignoring  the 
proportionality  factor,  Ppw  =  is  plotted  in  Figure  2  (small  circles)  for  a  range  of  values  of 
0^2. 


2.2  Energy  in  the  scattered  field 

No  loss  to  absorption  is  yet  included  in  the  modeltt.  (When  loss  is  included  (§3.6),  the 
conclusions  are  found  to  be  unaltered.)  Consequently  the  energy  in  the  scattered  field  Pscat  =  1- 
A^,  which  is  plotted  (as  stars)  in  Figure  2,  is  the  difference  between  the  incident  (unit)  power 

and  Ppw.  Lastly,  the  ratio  Pscat /Ppw.  is  plotted  as  plus  signs.  The  fields  have  equal  energies  at 
=  0.7  and  the  ratio  rises  to  1.7  at  c^=  1.  Large  values  of  severely  degrade  the  quality  of 

ultrasound  images.  From  the  theory  of  radiation  pattern  distortion  caused  by  random  phase 
errors*'  **  o|  must  fall  in  the  neighborhood  of  lO'^  or  less  for  reasonably  acceptable  beam 

quality.  Thus  practical  interest  lies  to  the  extreme  left  of  Figure  2,  which  implies  a  very  weak 
scattering  process  in  which  most  of  the  energy  passes  through  each  screen  unscattered. 


23  Angular  width  of  the  scattered  field 

The  scattered  field  is  the  difference  between  (2)  and  (3): 


fscat(u)  *  Jp(<t»  -  E<|))  -  (<j>2 .  E4>2)/2  -  j((|>3  .  E4>3)/3!+  -  E|)4)/4!  -)]  cj>““  dx 


[2n»l  -  E(j>n)  M!]  eikx»  dx 

[1  -  A  +  Sn.!  eHau  dx 


(5) 


fscat(u)  is  a  function  of  the  random  process  <►(*)  and  therefore  is  also  a  random  process.  Its 
power  Pscat  =  l-A^.  Its  angular  width  is  found  from  its  angular  power  density  spectrum  S(u). 
S(u)  is  proportional  to  the  expectation  of  the  square  magnitude  of  the  scattered  field,  E|fscat(u)i^. 
which  in  turn  can  be  obtained  from  the  Fourier  transform  F{»}  of  the  correlation  function 
calculated  from  the  measured  radiation  field,  i.e.,  S(u)  =  F{R(x)}. 

Wavefront  correlation  experiments  performed  in  in  vivo  breast^^  have  measured  the 
normalized  correlation  function  p(x)  =  R(x)/R(0).  From  these  data  the  correlation  function  is 
found  to  approximate  a  decaying  exponential  (Rgure  3),  Hence  we  model 


Absorption  is  added  later  in  §3.6,  where  it  is  found  that  its  inclusion  does  not  materially  alter  the  conclusions 
of  §2.2. 

*  Average  sidelobe  powerlevel  relative  to  main  lobe  =  ASL  =  where  Ng  is  the  no.  of  anay  elements^* . 
Specifying  -40  dB  for  the  required  ASL,  and  assuming  the  order  of  100  array  elements,  ^lO'^.  This  formula 
has  been  checked  repeatedly  in  adaptive  beamforming  experiments  in  radar  beginning  in  the  early  1980$'^. 
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p(x)  =  exp(-x/xo) 


(6) 


and  call  xq  the  e-^  correlation  distance#*,  xq  =  1.8,  3.3  and  3.3  mm  for  average  correlation 
curves  of  three  groups  of  asymptomatic  women:  premenopausal  dense  breasts,  premenopausal 
fatty  and  postmenopausal.  The  number  of  correlation  functions  averaged  were  35,  52  and  76. 
Other  reponed  measurements  are  in  the  same  neighborhood®*^®.  The  Figure  3  curves  and  the 
body  of  data  from  which  they  come  ^2  ^^re  obtained  over  long,  12-cm  paths,  which  corresponds 
to  many  correlation  distances  or  several  phase  screens  (§3.1).  Hence  it  is  possible  that  the 
correlation  evidence  leading  to  (6)  is  not  representative  of  a  single  screen.  This  is  checked  in 
§3.5  where  it  is  shown  that  the  angular  spectral  width  estimated  from  these  data  and  (6)  is 
representative  of  a  single  screen.  The  result  indicates  that  single  scattering  dominates,  which 
means  that  the  total  scattering  spectrum  observed  out  of  each  of  the  phase  screens  approximates 
a  first-order  spectrum.  This  further  implies  that  correlation  measurements  made  from  any 
thickness  of  tissue  disclose  xq. 

Based  on  (6)  the  unnormaliard  correlation  function  and  the  angular  power  spectrum  are 

R(x)-a|exp(-x/xo)  (7) 


S(u) 


1  +  (kuxo)2 


(8) 


and  the  half-power  width  of  the  scatter  power  spectrum  is 

Au  »  2/kxo  »  X/mio  (9) 

which,  for  a  frequency  of  3  MHz  (K  *  0.5  mm)  and  for  the  values  of  xq  given  above  for  the 
mean  correlation  functions  of  premenopausal-dense  and  premenopausal-fatty/postmenopausal 
breasts,  is  0.089  and  0.048  rad  or  2.7  and  5.1  deg.  Allowing  about  one  standard  deviation 
additional  spreads  about  the  means  the  predicted  halo  width  Au  increases  to  2-6  deg  at  3  MHz 
and  is  inversely  proponional  to  frequency. 

2.4  Estimate  from  breast  samples 

An  independent  prediction  of  halo  width  can  be  made  from  the  size  of  subcutaneous  fat 
globules  in  breast  tissue.  Fat  globules  have  random  sizes,  orientations  and  radii  of  curvature. 


As  distinct  from  the  FWHM  measure  of  correlation  distance.  The  relation  between  the  two  for  the 
exponential  correlation  model  (6)  is  =  1.38xo. 
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The  scale  is  predominantly  in  the  5  to  15  mm  range.  Each  globule  is  a  likely  candidate  for 
scattering  of  ultrasound  waves,  as  seen  in  a  radiologist's  schematic  of  a  breast  section  (Rgure 
4)13 

Referring  back  to  Figure  lb,  the  distance  between  scatterers  or  density  hills  shown  there 
is  nominally  a.  Scattered  energy  from  subcutaneous  fat  should  lie  largely  in  the  Fourio-  angular 
interval  Aup  -  Va,  5<  a  <15  mm,  which  at  3  MHz  is  0.033-0.1  rad  or  1.9  to  5.7  deg.  The  fact 
that  this  interval  is  substantially  identical  to  the  calculated  halo  width  of  2-6  deg  from 
correlation  data  (§2.3)  suggests  that  subcutaneous  fat  is  one  of  the  contributors  to  the  scattered 
field. 

2.5  Exponential  coefficient  vs.  frequency 

Time  delay  variations  should  be  independent  of  frequency  because  they  result  directly 
from  local  tissue  density  variations.  Frequency  dependence  manifests  itself  in  the  phase 
variance  where  o|  =  a)2oj2_  Consequently,  the  phase  variance  should  grow  as  the  square  of  the 

frequency.  Theory  also  predicts  that  xq  and  the  scatterer  size  a  are  related  by  a-  tcxq  by 
equating  the  two  calculations  of  scatter  spectrum  width,  Au  ■  X/jcxq  (9)  and  Aup  *  X/a  (§2.4). 
But  a  is  independent  of  frequency.  Therefore  the  correlation  function  should  be  independent  of 
frequency,  which  agrees  with  the  3  and  4  MHz  comparison  shown  in  Figure  3. 

2.6  Point  source  image 

When  the  composite  field  f(u)  =  fpw(Q)  +  fscat(u)  is  imaged,  whether  by  ^herical  lens  or 
transducer  (Figure  5),  the  image  (in  one  dimension)  in  the  focal  plane  contains  the  point  source 
and  the  apparent  source  distribution  of  the  scattered  field,  each  convolved  with  the  diffiraction 
pattern  of  the  lens  or  transducer.  The  target  is  reduced  in  amplitude  to  A  and  it  images  as  the 
diffraction  pattern  of  the  transducer,  with  half-power  width  X/L  The  angular  clutter  spectrum  is 
depicted  in  the  figure  as  a  continuum,  which  is  not  a  realistic  portrayal  because  it  is  a  Fourier 
transform  of  the  random  phase-perturbation  process  <|>(x);  hence  it  too  is  a  random  process  with 
peaks  and  valleys.  No  instantaneous  sample  function  is  so  smooth.  For  simplicity  the  clutter 
amplitude  spectrum  is  drawn  as  a  smooth  function,  which  may  be  interpreted  as  the  sq.  rt  of  the 
angular  clutter  power  spectrum  S(u).  The  half-power  width  is  X/rtxo.  The  clutter  power  is  l-A^. 
After  convolution  with  the  transducer  diffraction  pattern  the  half-power  width  of  the  clutter 
spectrum  equals  the  rms  sum  of  XItixq  and  X/L: 

Aus  =  [X/7Cxo)2  +  (X/L)2]1/2  =  X/Kxo®  X/L  (10) 

Figure  6  pictures  other  images.  Pan  a  is  the  image  without  the  phase  screen.  The  unit- 
amplitude  plane-wave  image  is  spread  to  the  diffraction  pattern.  It  consists  of  a  main  lobe  of 
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width  VL  plus  sidelobes.  L  is  the  size  of  the  aperture.  The  image  with  the  screen  is  part  b;  it  is 
the  same  as  the  image  in  Figure  5.  Its  width  is  •  TJkxq  for  high  resolution  imaging,  for  then 
L»7Cxo.  It  is  also  possible  for  the  opposite  condition  to  hold,  i.e.,  tdcq  >  L,  which  occurs  when 
the  nominal  scatterer  size  exceeds  the  aperture  size  of  the  imaging  system.  In  this  case  the 
"scatterer"  is  a  "large"  surface  (relative  to  the  transducer)  such  as  a  wedge  which  refracts  the 
incident  wave  and  bends  the  beam.  The  image  is  then  displaced  from  its  original  orientation  (in 
the  u-axis)  by  the  refraction  angle  %.  Part  c  depicts  the  case  of  wedge  refraction  plus  weak 
scattering.  A  realistic  "wedge"  would  have  rounded  rather  than  planar  facets  which  would 
change  the  focal  length  of  the  beam  and  therefore  defocus  the  image  (d).  For  the  pure  scattering 
problem  addressed  in  this  paper  we  ignore  refraction  and  select  Aus  *  X/tcxq. 

2.7  Resolution 

It  is  evident  that  lateral  resolution  is  unaffected  by  the  scatter  field;  it  remains  the  width 
X/L  of  the  mainlobe  of  the  diffraction  pattern  of  the  aperture  unless  the  scattered  energy 
obscures  the  main  image  lobe^^.  When  this  happens  lateral  resolution  and  contrast  vanish. 
This  phenomenon  is  a  direct  consequence  of  increasing  the  path  length  (§3.4  and  eq.(13)) . 

2.8  Target  to  clutter  power  ratio 

Figure  6  is  an  angular  amplitude  spectrum.  While  the  scatter  spectrum  is  pictured  as  a 
smooth  continuum,  any  sample  spectrum  is  itself  a  random  function  with  peaks  and  valleys 
because  the  spectrum  is  the  Fourier  transform  of  the  random  phase-screen  process.  The  target 
can  be  discerned  so  long  as  it  well  exceeds  the  peaks  of  the  scatter  spectrum.  The  measure  of 
target  detectability  is  the  target  contrast  or  target  to  clutter  power  ratio  TCR,  which  can  be 
estimated  for  a  square  2-D  array  of  size  L  as  follows. 

The  target  and  the  clutter  fields  spread  in  two  dimensions  as  they  propagate  (Hgure  7a). 
Approximate  the  angular  power  spectrum  by  the  sum  of  two  right  circular  cylinders  (part  b), 
pictured  in  cross  section  as  rectangles  in  part  c.  The  first  cylinder  is  the  target  spectrum.  Its 
power  A2  spreads  over  a  circle  of  diameter  X/L.  Its  height,  therefore,  is  (AL/X)2(4/7i).  The 
second  spectrum  is  the  clutter  region  or  halo  S(u)  induced  by  the  target.  Its  area  is  l-A^  and  its 
width  is  X/TtxQ.  Its  height  is  (l-A2)(jcxo/X)2(4/jt).  Hence  the  target  to  clutter  power  ratio  within 
the  halo  is 


TCR  =  ^  ^  .jj. 

(l-A2)(7tX0)2  O^(7ixo)^  C02o|(71X0)2 

where  the  approximation  =i>  A  =  1  and  1-A-  =  (^.  Eq.(l  1)  shows  that  TCR  for  a  single,  thin 
screen  is  quadratic  with  2-D  transducer  size  (linear  with  area)  and  inversely  proportional  to  the 
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time-delay  variance  of  the  screen  and  the  square  of  the  frequency  and  the  correlation  distance. 
The  frequency  dependence  of  TCR  enters  only  through  because  the  density  and 

compressibility  variations  in  the  tissue  induce  time  delay  variations  that  are  independent  of 
frequency  (§2.5). 

2.9  Target  dynamic  range 

From  TCR  (11)  we  can  determine  the  smallest  target  that  can  be  detected  in  the 
interference  background  and  hence  the  target  dynamic  range  (TDR),  which  is  the  ratio  of  the 
largest  to  the  smallest  detectable  target.  Eq.  (11)  is  the  ratio  of  the  target  intensity  to  the 
average  intensity  of  the  clutter  field  within  the  halo.  Because  the  field  develops  as  a 
bandlimited,  noise-like  random  process,  upward  deviations  of  10  dB  from  the  mean  are  not 
unlikely  (order  of  1%  occurrence  frequency)!^.  Hence  other  targets  within  the  halo  that  are  10 
dB  greater  than  the  mean  clutter  intensity  will  not  be  discernible  amidst  the  high  peaks  in  the 
clutter  spectrum.  At  least  a  few  more  dB  are  required  for  unambiguous  single-look  target 
detection,  although  with  multiple  looks  10  dB  is  sufficient.  Thus  suitable  contrast  between  the 
target  image  and  the  clutter  background  within  the  halo  requires  a  minimum  TCR  of  about  10 
dB.  In  other  words,  TDR  within  the  halo  is  about  an  order  of  magnitude  smaller  than  TCR. 

3.  MULTIPLE  ADJACENT  SCREENS 

The  input  to  the  first  screen  is  the  unit  amplitude  plane  wave  shown  in  Figure  1.  The 
emergent  plane  wave  plus  the  scattered  field  (5)  is  the  input  to  the  second  screen,  the  output  of 
which  is  the  input  to  the  third  screen,  etc.,  as  shown  in  Hgure  8.  The  emergent  radiation  is 
pictured  there  as  angular  (Fourier)  spectra  with  the  u-axis  parallel  to  each  screen. 

3.1  Screen  thickness  and  number  of  screens 

The  number  of  screens  N  is  the  propagation  depth  D  of  a  target  divided  by  the  screen 
thickness.  The  infinitely-thin  screen  model  calls  for  as  small  a  thickness  as  other  considerations 
allow.  Because  of  the  assumed  independence  of  the  phase  distribution  from  screen  to  screen, 
screen  thickness  can  be  no  smaller  than  scatterer  size  a.  This  is  also  reasonable  on  physical 
grounds,  for  otherwise  the  scale  of  the  scattering  stracture  would  be  altered  by  choice  of  screen 
dimension.  For  these  reasons  the  nominal  scatterer  size  is  taken  as  the  lower  bound  on  screen 
thickness.  Thus,  for  a  propagation  depth  D  of  6  cm  and  a  nominal  scatter  size  a  =  1  cm,  the 
number  of  screens  N  *  6. 

3.2  Output  angular  spectrum 
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The  clutter  build-up  is  illustrated  in  Figure  8.  Its  angular  spectrum  is  derived  as  the 
output  of  a  cascade  of  N  linear  filters,  each  with  angular  impulse  response  hn(u),  n  =  1,2...N. 
The  first  screen  has  impulse  response  hi(u)  =  AiS(u)  +  si(u),  which  is  the  total  spectrum  out  of 
the  Hrst  screen.  Ai5(u)  is  a  Dirac  5-funcdon  of  amplitude  Ai  (power  A^)  located  at  the  origin 
in  the  angular  spectrum.  si(u)  is  a  continuous  scatter  spectrum  of  power  1-Ai^»  a|,  the  phase 
variance  of  the  screen.  The  total  power  is  unity,  as  it  must  be,  because  no  losses  have  yet  been 
assumed.  The  width  of  si(u)  is  Au  ■  X/Jtxo  (9). 

Screen  2  has  impulse  response  h2(u)  s  A25(u)  +  S2(u).  Its  angular  spectrum  is 
hi(u)*h2(u),  where  ♦  denotes  convolution.  It  consists  of  three  parts;  a  plane  wave  of  amplitude 
A1A2,  the  sum  of  two  primary  scatter  spectra,  A2Si(u)  +  AiS2(u),  and  a  secondary  scatter  field 
si(u)  •  S2(u).  The  first  primary  scatter  spectrum  arises  from  the  si(u)  passing  unchanged  in 
shape  through  the  screen  but  reduced  in  amplitude  by  A2,  the  plane-wave  voltage  gain  of  screen 
2;  the  second  is  the  scattering  AiS2(u)  of  the  plane  wave  incident  upon  screen  2.  The  power  in 
the  primary  scatter  spectrum  is  A22(l-Ai2)  +  The  power  in  the  second  order 

scatter  spectrum  is  (1-Ai2)(l-A22).  The  total  power  in  the  emergent  field  is  again  unity 
(Ai2A22  +  A22[1-Ai2]  +  +  [1-Ai2l[l-A22]  »  1). 

The  angular  spectrum  fix>m  screen  3  it  is  hi(u)*h2(u)*h3(u).  It  consists  of  a  plane  wave 
AiA2A35(u),  a  primary  scatter  field  A3[A2Si(u)  +  AiS2(u)]  +  AiA2S3(u),  a  secondary  field 
[A2Si(u)  +  AiS2(u)]  •  S3(u)  +  A3[si(u)  •  S2(u)]  and  a  tertiary  scatter  field  si(u)  •  S2(u)  •  S3(u). 
The  primary  scatter  power  is  A32A22(1-Ai2)  +  A32Ai2(1-A22)  +  Ai2A22(1-A32).  The 
secondary  scatter  power  is  A22(1-Ai2)(1-A32)  +  Ai2(l.A22)(l-A32)  +  A32(1.A22)(1-Ai2)  and 
the  tertiary  scatter  power  is  (l-Ai2)(l.A22)(l-A32). 

After  N  screens  the  angular  spectrum  is  the  N-fold  convolution  hi(u)*h2(u)**—*hN(u). 
It  contains  a  5-function  representing  the  emergent  plane  wave.  Its  amplitude  is  AiA2'***Ai<f 
which  we  will  call  AN,  where  A  is  the  geometric  mean  of  the  Aj.  This  is  the  target  All  other 
components  of  the  output  spectrum  add  to  create  the  signal-induced  clutter  spectrum. 

The  power  spectrum  of  the  output  primary  clutter  field  in  Figure  8  is  the  sum  of  the 
primary  power  spectra  and  similarly  for  the  higher  order  scattered  fields.  Fortunately  this 
complicated  picture  can  be  simplified  considerably. 

3.3  Output  power  spectra 

Schema  of  cross  sections  of  the  target  spectrum  and  the  first  three  scatter  spectra  are 
illustrated  in  Figure  9.  The  output  target  power  Pt  =  Ai  ]2  =  a2N.  The  target  power  is 

spread  in  both  dimensions  over  the  diffraction  pattern  of  the  transducer.  The  bulk  of  the  power 
is  in  the  main  lobe.  Its  beamwidth  is  X/L  and  its  angular  cross  section  is  (7i/4)(VL)2-  Within  the 
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beam  the  power  density  is 


4A2N 

Jt(X/L)2' 


Its  spectrum  is  shown  schematically  in  the  upper  part  of 


Figure  9. 

The  Si(u)  are  independent  spectra  because  they  are  derived  from  independent  sample 
functions  of  the  random  phase-delay  process  defining  the  inhomogeneous  medium.  They  are 
continuous,  noise-like  terms  with  zero  means  and  their  powers  add.  The  same  is  true  for 
convolutional  terms  such  as  Si(u)  •  sj(u),  Si(u)  *  sj(u)  *  Sk(u),  etc.  What  is  peninent  for  each 
such  term  is  its  width  and  its  power  and  to  a  much  lesser  extent  its  shape.  All  Si(u)  have  the 
same  nominal  width,  power  and  shape.  Therefore  we  identify  all  SiCu)  with  a  generic  first  order 
random  scatter  spectrum  s(u),  all  Si(u)  •  sj(u)  with  a  generic  second  order  scatter  spectrum  s(u)  • 
s(u),  and  so  forth.  The  half-power  width  of  s(u)  is  X/jtxo.  Variances  add  under  convolution. 
Therefore  the  width  of  s(u)  •  s(u)  is  V2X/Jcxo.  The  widths  of  higher  order  convolutions  are 
VsX/rcxo,  V4X/71X0,  etc.  Areas  multiply  under  convolution.  The  power  in  s(u)  is  l-A^,  in  s(u)  • 
s(u)  it  is  (1-A^)^,  etc.  One  s(u)  term  arises  from  the  first  screen,  two  s(u)  terms  weighted  by  A 
from  the  second,  three  A^s(u)  terms  from  the  third,  and  so  forth.  Finally,  N  s(u)  terms  emerge 
from  the  final  screen  with  a  weight  of  A^**.  The  power  in  each  term  is  (1-A2)A2(N-1),  Because 
the  N  s(u)  terms  are  independent,  their  powers  add.  Therefore  the  power  in  the  first  order 
scatter  spectrum  Pci  =  N(1-A2)A2(N"1). 

Similarly  the  power  in  the  second  order  spectrum  Pc2  =  j(l-A2)2A2(N-2),  PC2«  Pci 

because  A  must  be  close  to  unity  to  achieve  adequate  TCR  for  useful  imaging.  For  the  third 
order  Pc3  *  ^^j(l-A2)3A2(N-3)«  and  so  on. 

The  power  density  is  estimated  in  the  same  way  as  was  done  for  the  target  The  power 
in  the  first  order  scatter  field  Pci  is  spread  over  the  scattering  angle  of  s(u),  which  is  X/itXQ, 
within  which  the  power  density  is  Pci/(^icx(j)2(j^4),  and  similarly  for  the  higher  order  scatter 
fields.  The  power  in  each  successive  scatter  order  decreases  by  the  factor  (l-A^).  Because  A 
is  close  to  unity,  this  ensures  that  the  first  order  scatter  field  dominates. 


3.4  Target  to  clutter  power  ratio 

Combining  the  terms  above,  the  target  to  clutter  power  ratio  for  N  screens  is 


TCR  = 


l^a^n 


(12) 
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£q.l2  can  be  simplified  considerably.  Figure  10  plots  TCR  vs.  »  l-A^  within  the  "clutter" 
region  for  L  =  25  mm,  Ttxo  =10  mm  and  N  =  6,  the  number  used  in  the  earlier  example.  The 
miniiniim  useful  TCR  must  be  larger  than  10  dB  because  the  target  dynamic  range  is  about  10 
dB  smaller  than  TCR  (§2.9).  Thus  the  maximum  phase  variance  per  screen  is  about  0.01  sq  rad 
for  unambiguous  detection.  Thus,  only  very  small  phase  variance  per  screen  can  be  tolerated 
for  otherwise  the  target  rapidly  gets  lost  in  the  background  as  grows.  Therefore  we  can 
assume  oj  »  10-2  and  «  1  -  Next,  Na  *  D  and  a  =  tcxq.  Also  <5^  =  Therefore 

(12)  becomes 


TCR 


ii«l 


n 


j 


ff \ 

-1 


7a,DcoW^ 


(13) 


Not  surprisingly  (13)  has  the  same  form  as  (11),  the  single  screen  TCR,  because  of  the 
dominance  of  the  primary  scatto*  field.  In  essence  both  are  weak  scattering  solutions. 

Eq.l3  predicts  that  TCR  is  inversely  proportional  to  propagation  distance.  9in  vitro 
TCR  measurements  were  made  at  the  Univ.  of  Rochester  laboratory  of  Prof.  Robert  Waag  from 
five  samples  of  different  thicknesses  of  excised,  frozen  and  defrosted  breast  tissue  used  in 
Ref.  10.  For  each  measurement,  a  tissue  sample  was  fixed  in  a  holder  within  a  water  tank. 
Propagation  was  one-way  from  a  pulsed  3.75  MHz  hemispherical  (virtual  point-)  source, 
through  the  water  and  tissue  path,  u>  a  1-D  linear  array  in  which  waveforms  were  recorded  at 
each  element  The  array  was  mounted  in  a  computer-controlled  2-D  posidoner  and  moved 
along  the  orthogonal  axis  so  as  to  form  a  92x46  sq  mm,  2-D  synthetic  array.  Waveforms 
measured  at  each  element  position  of  the  synthetic  2-D  array  were  imaged.  Figure  1 1  shows  a 
1-D  maximum-projection  cut  through  such  an  image.  TCR  measurements  were  made  from  such 
cuts.  The  central  lobe  is  the  image  of  the  source.  Clustered  around  it  are  discrete  multipath 
arrivals  from  the  source.  They  image  as  false  targets.  The  broad  clutter  background  is  the 
isotropic  scatter  field.  TCR  is  the  squared  reciprocal  of  the  estimated  clutter  level  at  the  origin. 
The  9  TCR  data  points  are  shown  on  the  log-log  plot  of  Figure  12a.  Each  point  is  the  average 
of  measurements  made  from  a  minimum  of  6  cuts.  The  straight  line  drawn  in  the  figure  has  the 
predicted  slope  of  -1.  With  the  exception  of  one  point  -  the  outlier  at  30  mm  depth  -  the  data 
approximate  the  predicted  behavior. 

3.5  Relations  between  single  and  multiple  screen  parameters 
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Correlarion  distance:  A  signifiqant  amount  of  multiple  scattering  is  precluded  because 
has  to  be  small  (-  lO^  sq  rad)  for  successful  imaging.  This  dominance  of  the  primary 

scatter  field  means  that  the  spectral  width  does  not  increase  with  propagation  path  length. 
Consequently,  xq  is  a  function  of  scatterer  size  only  and  not  the  thickness  of  the  propagation 
medium  in  which  it  is  measured.  Experimental  data  support  this  conclusion.  Measured  halo 
widths  from  different  thicknesses  of  a  single  in  vitro  breast  sample  used  in  Ref.  10  are  plotted  in 
Figure  12b.  The  angular  width  is  constant  over  a  thickness  range  from  15  to  40  mm.  Thus 
correlation  measurements  made  in  a  given  type  of  tissue  at  one  depth  are  pertinent  for  any  other 
depth. 

Phase  variance:  in  (13)  is  the  phase  variance  of  a  single  screen.  A  measurement  in 

tissue,  however,  is  never  across  a  single  screen  and  therefore  the  relation  of  measured  phase 
variance  data  to  must  be  determined.  We  can  view  the  sequence  of  screens  in  the  multi¬ 
screen  model  as  a  cascade  of  linear  filters  each  with  phase  delay  <|>i(x)  or  time  delay  ti(x).  The 
phase  (and  time)  delays  add  across  a  group  of  such  circuits.  Thus  a  phase  measurement  across  N 
screens  yields  <l>i(x)-H|t2(x)-H}>3(x)+...  =  (|)j(x)  and  because  of  the  assumed  independence  and 

common  statistics  of  the  <|>i(x)  the  variance  is  which  accounts  for  the  factor  of  N  in  (13). 

3.6  Absorption 

To  incorporate  absorption  in  the  model,  an  exponential  loss  exp(-(xz)  is  appended  such 
that  the  total  field  through  each  screen  of  thickness  a  is  attenuated  by  the  factor  exp(-aa)  where 
a  is  assumed  independent  of  positiion.  The  gain  of  the  ith  screen  to  the  plane-wave  component 
of  the  field  changes  from  Ai  to  a1  =  Aiexp(-(Xfl).  Output  target  power  Pj  =  [H^l  Aj^  =  A^^ 
(§3.3)  changes  to  aJ  ]2 
=  A2Nexp(-2Noui). 

The  clutter  is  modified  similarly.  All  scatter  field  terms  in  §32,  calculated  for  a  lossless 
model,  are  attenuated  by  the  factor  exp(-an).  The  first  order  spectrum  Pci  =  N(1-A2)A2(N*1) 
changes  to  N(l-A2)e-2aflA2(N-l)e*2(N.l)aa  =  N(l-A2)A2(N-l)e*2Nao.  Target  to  clutter  ratio 
remains  unchanged  for  the  primary  scatter  spectrum,  and  similarly  for  the  higher  order  spectra. 

The  same  is  true  for  the  ratio  Pjcat  /PpwJ  absorption  does  not  affect  the  balance  between 
the  scattered  and  plane  wave  components.  Further,  the  argument  holds  at  every  firequency  and 
therefore  frequency  selective  absorption,  which  alters  the  spectrum  of  a  deeply  propagating 
wideband  pulse,  does  not  affect  the  ratio  of  the  scattered  to  the  plane  wave  energy. 

Also  unaffected  by  absorption  are  the  relative  strengths  of  the  components  of  the 
scattered  field  that  arise  from  different  depths.  This  is  because  all  components,  whether  plane 
wave  or  scattered  field,  are  attenuated  identically  by  the  time  they  reach  the  skin.  To  illustrate 
this  phenomenon,  consider  primary  scattering  arising  from  the  nth  screen,  which  means  at  a 
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depth  of  na.  The  plane  wave  component  that  initiates  this  field  is  attenuated  by  absorption,  as  it 
travels  through  n-1  screens  toward  the  skin  surface,  by  the  factor  A2(n-l)e-2afl(n-l).  the  nth 
screen  the  scattered  field  s(u)  is  attenuated  by  the  factor  0*20^.  as  this  field  propagates  to  the 
skin  surface  it  suffers  additional  attenuation  in  the  remaining  N-n  layers  by  the  factor  A2(N-n)e. 
2(N-n)afl,  The  emergent  scattered  field  is  the  product  {A2(n-l)e*2afl(n-l))  {s(u)e*2ao)  {A2(N.n)e* 
2(N-n)«)  =  A2("-I*N.n)e-2aa(n.I+UN.n)5(„)  =  A2(N-l)e-2N<.o,{„),  j„ 

other  words  all  scattered  fields  arising  from  the  same  target  reach  the  surface  with  the  same 
strength.  Consequently  the  scattered  field  measured  at  the  surface  is  a  composite  of  the 
scattering  all  along  the  path  fiom  target  to  surface. 

3.7  Parameters  obtained  from  measurements 

Finally,  we  rewrite  (13)  in  terms  of  measurable  parameters.  Assume  that  measurements 
of  correlation  distance  and  time  delay  variance  are  made  across  a  body  of  length  »  a, 
where  a  is  the  nominal  thickness  of  a  phase  screen.  Call  the  measurements  Xom  and  The 
parameters  needed  for  analysis  arc  the  single-screen  values  Xq  and  or  O?.  Both  can  be 
obtained  from  the  measured  data.  As  discussed  above  (§3.5)  Xo  =  Xom  because  of  the 
dominance  of  primary  scattering.  The  single  screen  phase  variance  a^2  calculated  from  a^2  - 
shown  in  §3  J,  and  0^2  -  ©Zof.  Nm  is  the  number  of  screens  ot  scatterer  sizes  a  in 
Dm.  -  flq|)m/Dm  is  the  single  screen  phase  variance.  Eq.(13)  can  then  be  rewritten  as  a 
function  of  the  known  parameters  of  the  imaging  system  (ca,  L),  the  propagation  depth  during 

imaging  (D),  the  measured  tissue  parameters  (xom.  or  and  the  propagation  depth 
during  measurement  (Dm): 


Na|(7txo)2  D(icxo)2(a4i/Dm)  “  Da)2(nxom)2(a4/Dm)  ^  ^ 

where  the  last  approximation  uses  a  -jcxom  and  a^i  =  The  right-hand  expression  is 

the  most  generally  useful  form. 

The  numerator  parameter  and  first  two  denominator  parameters  are  design  parameters. 
The  next  two  denominator  parameters  arc  measured  properties  of  the  tissue.  The  last  is  the 
length  of  the  tissue  sample  in  which  tissue  parameter  measurements  are  made.  Thus,  time  delay 
variance  and  correlation  distance  arc  the  primary  tissue  measurements  required  for  prediction  of 
TCR  in  a  purely  scattering  medium. 


3.8  Limit  to  resolution 

In  §2.7,  the  lateral  resolution  was  found  to  remain  X/L,  the  width  of  the  diffraction 
pattern,  and  unaffected  by  the  scatter  field  so  long  as  the  scattered  energy  does  not  obscure  the 
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main  image  lobe.  When  detection  vanishes  so  do  lateral  and  contrast  resolution,  as  was 
observed  in  12-cm  path-length  observations  reported  in  Refs.7  and  8. 

4.  OTHER  RELATED  TOPICS 

4.1  Near  fleld 

The  theoretical  development  of  the  scatter  field  used  a  plane  wave  or  far-fleld  model 
whereas  ultrasonic  imaging  instruments  are  focused  in  the  near  field.  There  are  two  separate 
matters  here.  First,  the  physics  of  scattering  has  nothing  to  do  with  the  focal  properties  of  the 
measuring  instrument  and  hence  is  unrelated  to  whether  or  not  the  instrument  is  focused. 
Second,  the  scatter  spectrum  from  which  spectral  width  was  calculated  extends  only  a  few 
beamwidths  from  the  focus  and  hence  is  satisfactorily  approximated  by  the  far-field  Fourier 
transform  calculation  of  the  scattered  field  in  (1).  Thus  the  plane  wave  model  is  acceptable. 

4.2  Two-way  propagation 

The  model  deals  with  both  point  source  radiation  propagating  one-way  to  a  receiving  array 
and  to  the  two-way  problem  of  an  isolated  point  target  when  illuminated  by  an  arbitrary 
transmittCT.  It  does  not  cover  the  arbitrary  target  distribution  with  arbitrary  insonification.  In 
principle  the  method  can  be  extended  to  any  source/target/transmitter  configuration  and  in  this 
paper  two  practical  and  easily  describable  cases  are  given. 

One  case  is  the  well-focused  transmitter.  Only  the  target  is  illuminated  with  full 
strength.  All  other  targets  receive  sidelobe  Ulumination.  Therefore  the  scatter  level  from  the 
target  is  higher  than  the  otho-  scattered  fields  by  the  sidelobe  level  of  the  transmitter.  Those 
fields  are  not  first  order  competitors  to  the  primary  scatter  background  from  the  target  and  in 
general  can  be  ignored. 

In  another  important  case  the  transmitter  focus  is  degraded  by  the  inhomogeneous 
medium  from  X/L  to  X/tcxq,  the  width  of  the  primary  clutter  spectrum,  where  L»Xo.  Now  all 
targets  within  a  2-D  circle  of  diameter  VitXQ  are  illuminated  nearly  equally  and  contribute 
equally  to  the  clutter  background.  For  the  sake  of  argument,  let  there  be  L/ttxo  targets  per 
diameter  (one  per  beamwidth)  of  nominally  the  same  size  (such  as  Rayleigh  scatterers)  so 
insonified.  Then  the  number  of  illuminated  scatterers  in  a  circle  of  diameter  L/itxo,  which  is 
(7i/4)(L/7:xo)2,  is  the  factor  by  which  the  clutter  power  increases.  The  total  received  clutter 
power  is  spread  in  two  dimensions  over  roughly  twice  the  area  (diameter  is  V2Vxxo  because  the 
received  spectrum,  being  the  convolution  of  the  scatter-source  spectrum  on  transmit  and  the 
medium-induced  spectrum  on  receive,  is  V2  times  larger  in  diameter).  The  result  is  a  change  in 
TCR  by  the  factor  (8/7t)(7ixo/L)^. 
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A  more  general  case  is  an  arbitrary  target  distribution  T(u)  illuminated  by  a  field  f(u) 
from  a  transmitter  radiating  through  an  inhomogeneous  medium.  Its  scatter  Held  is 
T(u)f(u)*s(u),  the  convolution  of  the  medium  response  to  the  target  reradiation.  Further,  let  the 
target  distribution  be  discrete  and  described  by  T(u)  =  ETk5(u-uk)  where  the  Tjc  are  the 
amplitudes  of  targets  located  at  u=uic.  The  scatter  field  is  s(u)*f(u)LTk5(u-Uk)  =  s(u)«I^kfkS(u* 
Uk)  =  2Tkfks(u-Uk)  where  fk  =  f(u-Uk). 

4.3  Limitations  of  the  model 

Infinite  statistically  homogeneous  phase  screen:  This  model  of  spatially  stationary 
statistics  is  satisfactory  for  weak  scattering  in  a  single,  large  tissue  bed.  It  fails  to  account  for 
interfaces  between  tissue  beds  across  which  a  sound  velocity  change  takes  place.  In  other 
words,  the  model  ignores  all  multiple  angular  arrivals  at  the  transducer  array  from  a  common 
source  which  could  arise  from  refraction  at  an  interface  or  reflection  from  a  strongly  reflecting 
surface.  Such  arrival  wavefronts  will  interfere  in  the  aperture  and  appear  to  the  imaging  system 
as  multiple  sources,  each  with  its  own  scattering  halo. 

The  model  is  expected  to  be  realistic  in  a  large,  reasonably  homogeneous  organ  such  as 
liver  and  for  the  purely  scattering  conqmnent  in  breast  from  subcutaneous  and  intramammary 
fat  that  fill  much  of  most  breasts.  High  density  tumors  and  liquid-filled  cysts  will  refraa 
wavefronts  and  introduce  image  artifacts.  It  is  hoped  that  these  will  submit  to  modeling  at  a 
later  time. 

Array  size:  The  model  is  appropriate  for  arrays  that  are  not  too  small  to  perceive  the 
spatial  fluctuations  in  <|>(x),  for  then  a  statistical  description  of  the  medium  is  inappropriate,  nor 
so  large  as  to  violate  spatial  stationarity  (fluctuation  statistics  no  longer  constant  across  the 
array).  The  smaller  limit  likely  occurs  when  L  «  a.  In  the  other  extreme  the  array  not  only 
sees  a  stochastic  field  of  random  spatial-densiQr  fluctuations  but  also  discrete  sub-bodies  (cysts, 
tumors)  with  different  refractive  indices.  In  this  case  the  model  given  in  this  paper  pertains 
only  to  the  purely  scattered  portion  of  the  field,  to  which  it  can  be  applied,  but  not  to  the 
coherent  multipath  from  the  discrete  bodies. 

5.  SUMMARY  AND  OBSERVATIONS 

Inhomogeneous  tissue  removes  energy  from  a  wave  both  by  incoherent  scattering  and 
coherent  refraction  and  reflection.  The  incoherent  scattered  field  is  studied  in  this  paper.  The 
modeling  technique  is  to  slice  the  tissue  body  into  thin,  random  phase  screens,  calculate  the 
properties  per  screen  and  combine  them.  Scattering  is  produced  by  random  speed  fluctuations 
in  tissue,  largely  due  to  density  and  compressibility  perturbations.  The  nominal  distance 
between  density  hills  is  identified  with  scatterer  size  a  =  ttxo,  where  xq  is  the  exponential 
correlation  distance.  Scanering  creates  a  halo  around  the  target  that  diminishes  the  image 
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contrast.  The  halo  width  is  *  TJa.  The  loss  in  contrast  is  measured  by  the  target  to  clutter 
power  ratio,  which  diminishes  with  the  correlation  distance,  the  propagation  length  D  in  the 
tissue  and  the  time-delay  variance  per  unit  distance,  and  the  square  of  the  frequency.  Lateral 
resolution  within  the  halo  is  unaffected  by  the  scatter  field;  it  remains  X/L.  L  is  the  array  size. 
However,  detection  of  weak  targets  within  the  halo  of  a  strong  target  requires  that  TCR  »1. 
When  TCR  is  small,  however,  detection  is  not  possible  and  the  preservation  of  the  lateral 
resolution  is  a  meaningless  vinue. 

The  obscuring  effect  of  the  halo  is  not  only  from  its  intensity  but  also  from  the  number 
of  lateral  resolution  cells  in  which  the  clutter  can  hide  targets.  This  number  is  the  ratio  of  the 
halo  cross  section  to  the  beam  cross  section,  which  is  (L/fl)2.  This  ratio  is  independent  of 
wavelength.  Scattering  occurs  at  all  depths  between  target  and  skin  surface.  The  principal 
scattering  mechanism  is  single  or  first-order  scattering  because  the  energy  removed  by  each 
scattering  hill  is  small. 

The  following  summarizes  comparisons  between  predictions  of  the  theory  and 
experimental  data. 

1.  Halo  width:  Based  on  the  size  of  fat  globules  the  theory  predicts  2-6  deg  at  3  MHz  (§2.4). 
Calculations  from  wavefiont  correlation  experiments  give  the  same  spectral  width  (§2.3). 

2.  e*^  correlation  distance  vs.  frequency:  The  experiment  shown  in  Figure  3  indicates 
independence,  which  is  what  the  theory  predicts  (§2.5). 

3.  TCR  vs.  propagation  depth:  Theory  predicts  inverse  linear  dq>endence.  Experiment  shown 
in  Figure  12a  indicates  approximately  the  same  behavior. 

4.  Halo  width  vs.  propagation  depth:  Theory  predicts  independence.  Experiments  show  the 
same  (Figure  12b). 
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Hguie  1.  (a)  Random,  stadsdcally  homogeneous  medium  modeled  as  sequence  of  adjacent 
random  phase  screens,  (b)  De^  of  single,  infinite  screen  showing  density  variadons 
(contours)  as  cause  of  local  perturbadons  of  dme  delay  through  screen.  Some  of  the  energy 
of  the  plane  wave  incident  upon  the  screen  emerges  as  a  plane  wave.  The  remainder  is 
scattered.  Shown  is  the  scattered  wavefront  from  one  scattering  center  in  the  screen. 
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Figure  2.  Emergent  energy  vs  phase  variance  -  (2pfst)2.  Single  phase  screen.  Ppw  ® 

a2  =  plane  wave  power,  Pjcat  *  =  scattered  power,  Pscat/Ppw  *  ratio  of  scattered  to 
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Figure  3.  (a)  Point  source  coirelation  functions  of  12-cm  thickness  in  vivo  breasts 
measured  2.5-3  cm  from  chest  wall.  Data  from  Ref  1.  (b)  Average  of  the  three 
curves  of  (a)  and  exponential  fit. 
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Ganglion. 


Figure  4.  Schematic  of  breast  section.  From  Ref.  13. 
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Figure  5.  Image  contains  diffraction  pattern  of  aperture  (image  of  emergent  plane  wave) 
and  angular  spectrum  of  the  scatter  field. 
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Figure  6.  Image  of  point  source  is  diffraction  pattern  of  imaging  apenure.  (a)  Without 
phase  screen,  (b)  Scatter  spectrum  caused  by  phase  screen.  Image  power  reduced,  (c) 
Linear-wedge  refraction  plus  weak  scattering.  Image  displaced  through  angle  Uf.  (d) 
Wedge  with  rounded  facets  defocuses  image. 
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Figure  7.  Angular  target  and  clutter  spectra,  (a)  Contour  plot  of  point  target  (center 
contours)  plus  widely  dispersed  scattered  field  and  ^screte  multipath.  The  coordinates  are 
azimuth  and  elevation.  In  vitro  breast  sample  (brs006)  is  4-cm  thick,  3.75  MHz.  92  mm  x 
46  mm  array.  Outer  contour  -15  dB,  Data  from  experiment  reponed  in  Ref.  10.  (b) 
Approximate  target  and  clutter  spectra  showing  heights  and  widths. 
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Figure  8.  Build-up  of  the  multiple-screen  angular  output  spectrum. 
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Figure  9.  Cross  sections  of  approximate  target  and  clutter  spectra  after  N 
screens.  2-D  array. 


Figure  10.  Ratio  of  point  target  intensity  Figure  11.  1-D  cut  at  >25  deg  through  single- 
in  image  to  clutter  background  induced  source  image  of  Figure  7a.  System  noise  level 

by  scattering,  is  the  phase  variance  •=  -35  dB.  Cut  is  through  region  with  minimum 
per  screen.  6  screens,  25  mm  array  size,  multiple  arrivals.  Central  lobe  is  target  image. 

10  mm  scatterer  size.  Plot  of  (13).  Adjacent  lobes  arc  multipath,  probably  refraction 

induced.  Low,  broad  background  is  isotropic 
scattering.  TCR  and  halo  width  are  estimated 
from  the  scatter  background. 
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Rgure  12.  (a)  TCR  and  (b)  halo  width  (deg)  vs.  thickness  of  in  vitro  breast  samples.  TCR 
shows  approximate  inverse  linear  dependence  (solid  line)  with  path  length.  Halo  width  shows 
independence. 


